Universal mechanism of spin relaxation in solids 
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We consider relaxation of a rigid spin cluster in an elastic medium in the presence of the magnetic 
field. Universal simple expression for spin-phonon matrix elements due to local rotations of the 
lattice is derived. The equivalence of the lattice frame and the laboratory frame approaches is 
established. For spin Hamiltonians with strong uniaxial anisotropy the field dependence of the 
transition rates due to rotations is analytically calculated and its universality is demonstrated. The 
role of time reversal symmetry in spin-phonon transitions has been elucidated. The theory provides 
lower bound on the decoherence of any spin-based solid-state qubit. 
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I. INTRODUCTION 

Understanding spin-lattice interactions has fundamen- 
tal importance for applications of magnetic phenomena, 
such as, e.g., magnetic relaxation, magnetic resonance, 
and decoherence of spin-based qubits. The studies of 
spin-lattice interactions are almost as old as the quantum 
theory of solids pLSiliSi Van Vleck* was among the first 
who attempted to compute the rates of phonon-induced 
transitions (in titanium and chromalum) from the first 
principles, that is, considering atomic wave functions in 
the crystal field and spin-orbit interaction. These cal- 
culations are very involved, and they are hardly pos- 
sible for more complicated materials such as magnetic 
molecules or clusters, each cluster containing several 
magnetic atoms and tens or hundreds of nonmagnetic 
atoms. 

On the other hand, for most materials the orbital mo- 
ment of electrons is quenched by a strong crystal field 
and it does not contribute to the magnetic properties. 
The latter are due to the spin S that is relatively weakly 
coupled to the orbital moment via the spin-orbit inter- 
action and thus it feels the crystal field indirectly. While 
microscopic calculation of the crystal-field Hamiltonian 
for S is difficult, one can use an alternative approach 5 6 
and start with a phenomenological expression containing 
all permitted by symmetry combinations of spin S. Some 
examples of the crystal-field Hamiltonian that results in 
the uniaxial, biaxial, and cubic magnetic anisotropy of 
the cluster are given below: 
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where e'"' with a = 1, 2, 3 are the unit vectors of the co- 
ordinate frame that is rigidly coupled with the symmetry 
axes of the cluster (in the following it will be called "lat- 
tice frame"). The advantage of this approach is that the 
form of Ha is uniquely determined by symmetry, while 
the phenomenological constants (D, E, C, etc.) can be 



easily obtained from experiment. The expressions for Ha 
above, as well as all other physically acceptable forms of 
Ha, possess a full rotational invariance, i.e., the form of 
Ha does not depend on the orientation of laboratory co- 
ordinate axes. The full rotational invariance implies con- 
servation of the total angular momentum (spin + lattice). 
Eqs. Q-© are also applicable to magnetic molecules or 
clusters, S being the total spin of the cluster formed by a 
strong exchange interaction between the individual mag- 
netic atoms. 

Lattice vibrations couple to the spin by modifying coef- 
ficients in Eqs. changing directions of the lattice 
vectors e^") , as well as by lowering the symmetry of the 
crystal field. Due to translational invariance, the phonon 
displacement field u enters the Hamiltonian only through 
its spatial derivatives dui / dx.j or their symmetric and an- 
tisymmetric combinations. Again, one can write down a 
phenomenological expression for the spin-phonon Hamil- 
tonian -ffs-ph that contains all terms permitted by the 
symmetry of a particular material. 7 However the general 
form of iis_ph usually contains too many different coef- 
ficients that are comparable with each other and cannot 
be measured independently. 

The problem of spin-lattice relaxation can be simpli- 
fied if one notices that longitudinal phonons have a larger 
sound velocity than the transverse phonons. Since the 
rate of one-phonon processes (emission and absorption 
of a phonon) is inversely proportional to the fifth power 
of the sound velocity, procecces involving longitudinal 
phonons can be neglected. The same is valid for mul- 
tiphonon processes such as the Raman process, because 
their rates contain even higher powers of the sound ve- 
locity. 

The terms of /f s _ p h due to transverse phonons can be 
split into two groups. The first group describes distor- 
sions of the lattice cell due to transverse phonons whereas 
the second group of terms describes local rotations of the 
lattice without distortion of the crystal environment of 
magnetic atoms. Whereas the first group of terms con- 
tains phenomenological coupling coefficients, the terms of 
i/s-ph due to local rotations are parameter free and are 
defined solely by the form of Ha- The significance of the 
latter was noticed in the past, 8 ^ 10 : 11 ! 12 : 13 ! 14 : 15 and dif- 
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ferent kinds of magnetoelastic problems have been con- 
sidered. However early applications of the theory have 
not included spin-lattice relaxation. Much later the re- 
laxation between the adjacent spin levels of the spin 
Hamiltonian of Eq. Q due to the parameter-free spin- 
lattice interaction arising from tilting of e^ 3 ) by trans- 
verse phonons was considered in Ref. 16. 

In general, processes due to the distortion of the lat- 
tice and those due to the local rotation of the lattice 
should result in comparable relaxation rates. Even in 
this case, the latter are of a fundamental importance be- 
cause they provide a parameter-free lower bound on the 
decoherence of any spin-based qubit. In the case of a 
magnetic molecule or a cluster in a solid, if the cluster 
is more rigid than its environment, it resists any distor- 
tions due to long-wave deformations of the solid. That is, 
S interacts only with the long-wave deformations of the 
crystal lattice that rotate the local frame (e^, e( 2 \ e( 3 )) 
as a whole, so that the spin-phonon interaction can be 
obtained from Ha without any phenomenological param- 
eters. The corresponding parameter-free description of 
spin-lattice relaxation becomes exact in this case. 

It has been noticediiiSii2iS£ that for an arbitrary spin 
Hamiltonian the parameter-free spin-phonon transition 
rates can be conveniently computed by switching to the 
lattice frame where the form of Ha is preserved, while the 
spin-phonon interaction is of a kinematic origin and it has 
a universal form that is independent from Ha- With this 
new method, one could easily calculate the relaxation 
rates between the tunnel-split states of Ha, while it was 
unclear how this problem could be solved by conventional 
methods in the laboratory frame. A striking feature of 
the lattice-frame approach is that both the calculation 
and the final result for the relaxation rate are universal 
and insensitive to the detailed form of Ha- One obtains 
the relaxation rate that is expressed via the quantities 
that can be directly measured in experiment, such as the 
tunnel splitting A. 

The aim of the present work is to investigate the rela- 
tion between the lattice-frame and laboratory-frame ap- 
proaches in more detail. We will show that also in the 
laboratory frame one can formulate a new method of cal- 
culating parameter-free spin-phonon rates that is similar 
to the lattice-frame approach and leads to the same re- 
sults. We will extend the theory by taking into account 
the magnetic field, including the case in which the tunnel 
splitting is solely due to the magnetic field. 

The body of this paper is organized as follows. In 
Sec. [n] we derive exact expressions for the spin-phonon 
interaction induced by rotations in the laboratory (Sec. 
Ill A|) and lattice (Sec. IIIB|I frames. In Sec. IIIII a simple 
universal formula for the spin-phonon matrix element is 
obtained for an arbitrary spin Hamiltonian in the pres- 
ence of the magnetic field. Equivalence of the laboratory 
and lattice frame treatments of the spin-phonon interac- 
tion is demonstrated. The advantage of our method is 
that it requires only the knowledge of the matrix elements 
of the operator S between the eigenstates of the Hamil- 
tonian, in contrast with the traditional method that re- 
quires the knowledge of a matrix element of a model- 



dependent product of the components of S. In Sec. IIVI 
spin matrix elements for transitions between spin states 
split by the crystal field or by the magnetic field are cal- 
culated analytically for Hamiltonians that arc dominated 
by the uniaxial anisotropy. We show that these matrix 
elements are parameter free and their field dependence is 
universal. At the end of this section we discuss the role 
of time-reversal symmetry. In Sec. we obtain univer- 
sal parameter-free formulas for the rates of spin-phonon 
transitions between tunnel-split states in the presence of 
the arbitrarily directed magnetic field. Cases of the tun- 
nel splitting due to the crystal field and due to trans- 
verse magnetic field is considered in Sees. IV Al and IV Bl 
respectively. Implications of our results for experiment 
are discussed in Sec. IVII Some illustrations for uniax- 
ial and biaxial spin models with the magnetic field, as 
well as in-depth study of the role of the time-reversal 
symmetry, are presented in Appendices. In particular, 
in Appendix [E] we show that the same expressions for 
the rates of relaxation between tunnel-split states can be 
obtained by traditional methods but the required effort 
greatly exceeds the effort of our new method. 

II. SPIN-LATTICE INTERACTION 

A. Laboratory frame 

In the absence of phonons, one can choose the coor- 
dinate system in which in Eqs. Q-© ejj = 8 a /3, i.e., 

Ha = —DS 2 etc. A transverse phonon, u(r), rotates 
the axes of the local crystal field, (e^\ e^- 2 \ ) . This 
rotation can be described by 

<50(r) =lv x u(r) (4) 

and it is performed by the (3 x 3) rotation matrix R, 

e (°0 _> R e ( Q ) , (5) 
i.e., ejg — ► R^'el, . For small 5<fi, one has 

R Q /3 = S a p - e Q /3 7 ^ 7 - (6) 

We now notice that due to the rotational invariance of 
Ha, the rotation of the local frame (e' 1 ), e' 2 \ e' 3 )) is 
equivalent to the rotation of the vector S in the oppo- 
site direction, S — » R _1 S. As is known^Si this rotation 
can be equivalently performed by the (25 + 1) x (25 + 1) 
matrix in the spin space, 

S -> RSR- 1 , R = e- lS . (7) 

The total Hamiltonian can be written in the form 

H = RH A R- l +H z + H ph , (8) 

where Ha is the crystal-field Hamiltonian in the absence 
of phonons, 

H z = -9A%H ■ S (9) 
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is the Zeeman Hamiltonian, and H p h is the Hamiltonian 
of harmonic phonons. 

In the above formulas, u and (f> must be treated as 
operators. Canonical quantization of phonons and Eq. 
® give 
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where M is the mass of the unit cell, N is the number of 
cells in the crystal, ekA are unit polarization vectors, A = 
ti,tQ,l denotes polarization, and lu^x = v\k is the phonon 
frequency. In application to rigid magnetic clusters, Eqs. 
HI U| l and Hllfl describe quantized long wave phonons in 
the elastic environment of the cluster. In the linear order 
in phonon amplitudes one obtains 



RHaR =HA+H s - P h, H s -ph — i 

The total Hamiltonian can be written as 
H = H + Hs-ph, 



H A ,S 



dep. 

(12) 

(13) 



where Hq is the Hamiltonian of non-interacting spin and 
phonons, 



Hn — Hq 



H, 



ph ) 



and 



H S = H A 
is the spin Hamiltonian. 



H 2 



(14) 
(15) 



B. Lattice frame 

The anisotropy Hamiltonian Ha is defined in the lat- 
tice frame and thus, in this frame, it is not changed by 
the lattice rotations. Unitary transformation to the lat- 
tice frame corresponds to the rotation of H of Eq. © by 
the angle 8(f): 



H' — R 'lilt — II \ + H' z + //.',,, 



Here 



and 



H'ph = R 1 HphR — H p h — i 



Hph, Sep 



(16) 
(17) 



H'z 



R ^HzR — Hz 



Hz, S 



8(f) 



Hz - 5Ms [Hx(50] -S, 



(18) 



where we have used [(A • S) , (B • S)] = iS- [A x B] . 

The full Hamiltonian in the lattice frame up to the first 
order in 5<f> is thus 



H = Hq + H 1 



s— ph J 



(19) 



where Hq is given by Eq. (|14|) and 

#s'- P h = - { H p h, 8(f) -S - gn B [Hx§4>] S. 
With account of the relation 



one obtains 



where 



S(P = -[H ph ,S(f>] 



H^-ph - ■ S, 



{l = 8if) + 7 [H x 8(f)] , 



(20) 



(21) 



(22) 



(23) 



and 7 = g[i B /h is the gyromagnetic ratio for S. 

Note that in the absence of the magnetic field the spin- 
lattice interaction in the lattice frame can be obtained 
by simply writing 17 i 18 i 19 i 20 H' = H - • S. The term 
— h8(f) • S is of kinematic origin: in the rotating coordi- 
nate frame the rotation is equivalent to the magnetic field 
H e fi=60/7 acting on the spin. The second term in Eq. 
(|2*3Tl describes the fact that the external magnetic field, 
which is constant in the laboratory frame, makes rotation 
in the lattice frame due to the transverse phonon. 



III. MATRIX ELEMENTS OF SPIN-LATTICE 
INTERACTION 

We study spin-phonon transitions between the eigen- 
states of Ho that are direct products of the spin and 
phonon states, 



(24) 



Here \ip±) are the eigenstates of H$ with eigenvalues E± 
(E + > E-) and \<f>±) are the eigenstates of H p h with 
energies E p h,±- For H s - p h linear in phonon amplitudes, 
the states \4> ± ) differ by one emitted or absorbed phonon 
with a wave vector k. Thus we will use the designations 

|0 + ) = |n k ), \(f>_) = \n^ + l). (25) 

Spin-phonon transitions conserve energy: 

E+ + £ ph ,+ = + £ p h,- (26) 

We calculate spin-phonon relaxation rates in both lattice 
frame and laboratory frame and show that the result is 
the same, as physically expected. 



A. Lattice frame 

We first calculate the matrix clement corresponding to 
the decay of the spin V+) ~* \^-) m l ne lattice frame. 
From Eq. I|22|l one obtains 



-^s-ph 



-M2_+-(^_|S|^ + >, (27) 
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where 



(28) 



To calculate the matrix element fi |_, it is convenient to 

step back and use the commutator form of S<fi given by 
Eq. (2TJ. One has 



iJph, 54> 



= (E ph ^-E pK+ )Sc/>_ + , (29) 



where 



= {<t>_\5<t>\<i> + ) . (30) 
It follows from the energy conservation, Eq. I|26(). that 

-Eph,- — -Eph,+ = E + — E- = hujQ. (31) 
Thus one finally obtains 

= iuj a 54>_ + + 7 [H x 5cf)_ + ] . (32) 

B. Laboratory frame 

To check the consistency of our method, let us now 
obtain the expression for the matrix element in the lab- 
oratory frame. Eq. Q12f) gives 



s— ph 



H A ,S 



(33) 



It is convenient to avoid explicitly working out the com- 
mutator in the spin matrix element. To this end, we add 
and subtract Hz- 



H A ,S 



i(ip- 



H Z ,S 



H S ,S 



(34) 



Now we can take into account that the states \ip±) are 

exact eigenstates of the spin Hamiltonian Hs of Eq. 1|15|) 
with energies E±: 



H S ,S 



V+) = - ^7+) <V»- |S] . (35) 



The Zeeman term in Eq. (|34|l can be done as in Eq. (|18fl . 
With the help of Eq. I|31|) one then obtains 



H 



s — ph 



*+) = -?>«-+ •(V-|S|V+>, (36) 



where f2 (.is given by Eq. I|32l) . We see that the spin- 

phonon matrix elements computed in the laboratory and 
lattice frames, Eqs. (|27() and i|36|) . are exactly the same. 



£„, K 



-20 



-40 



-60 



-80 



j i i i i i i i i i i i i i i i_i i i i_ 



m+1 /-. 
m— 1 J-* 



*\ m'—l 

»\ m = —m—k 

»N. m'+l 



Energy levels 
of the Mni 2 Ac cluster 

H = - DS 7 2 - H 7 S 7 



HJD = k = \ (first resonance) 



~i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — >- 1 — 1 — i — 1 — r 
10 -8 -6 -4 -2 2 4 6 8 10 

m 



FIG. 1: Spin energy levels of a Mn^Ac molecule for H x — 
and H z = D corresponding to the first resonance, k = 1. 



over conventional methods. Regardless of the explicit 
form of the spin Hamiltonian, only the matrix elements 
of the operator S need to be computed. Especially inter- 
esting is the case of strong uniaxial anisotropy, in which 
Hs nearly commutes with S z , so that the energy lev- 
els can be approximately described with the help of the 
quantum number m: 



E m = #m ) - gi^B H zm, 



(37) 



where Em is the contribution of the crystal field that 



satisfies E_ 



(A) 



Em^ ■ The structure of the energy levels 



for this model is shown in Fig. ^ The two levels to and 
to' are in resonance for the values of the magnetic field 



rr( res ) 
9^B n z,mm> 



to — to 



(A) 



m — m 



(38) 



For Mni2 the dominant term in the uniaxial anisotropy 
energy is E^' = -Dm 2 and thus gfJ- B H z ,mL' = ~D(m+ 
m') — Dk 7 where k is the resonance number (see Fig.^l. 
The level bias is given by 

W ee E m - E m , = g» B (H z - H<$J) (m' - m) 

= (to' - to) g\i B H z + hsJ^,. (39) 

Choosing dW/dH z > for certainty, makes m! > m in 
all subsequent calculations. 



IV. SPIN MATRIX ELEMENTS FOR 
TUNNEL-SPLIT STATES 



A. The two-state model: \S z \il> + ) 



The method of the computation of transition rates out- 
lined in the previous sections has significant advantage 



Due to the terms in Hs that do not commute with 
S z , the true eigenstates of Hs are expansions over the 
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complete |m) basis: 

s 



(40) 



-s 



If one neglects tunneling (i. e., hybridization of the states 
\m) and \m') on different sides of the potential barrier, 
see Fig. QJ, each pure state \m) should be replaced by 
|ra): 



m" — — S 



(41) 



where 



1 and all other coefficients are small. Hy- 



bridization of the states 



and 



can be taken into 



account in the framework of the two-state model 



(fhi 
(rh 



E„ 



-Ae llp , 



(42) 



where A is the tunnel splitting of the levels to and to' 
that can be calculated from the exact spin Hamiltonian 
Hs or determined experimentally and ip is a phase. If Hs 
expressed in terms of S z and S± is real, then the only two 
possibilities are if = 0, tt, so that e l(p = ±1. Diagonalizing 
this 2x2 matrix yields the eigenvalues 



E± = \ (E m 4 
The energy difference is 



E m , ± x/W 2 + A 2 



E+-E- = huj a = ^W 2 + A 2 . 



(43) 



(44) 



The corresponding eigenvectors can be represented in the 
form 



\1>±) = -i (C±e iv/2 \m) ± C T e~ lv/2 \fh') 
v2 ^ 



where 



C± = Jl± 



w 



Vw 2 



(45) 



(46) 



Far from the resonance, \W\ 3> A, the eigenstates and 
energy eigenvalues reduce to those of |to) and \fh') states. 
For W > and W > A one has \ip + ) = e lv/2 \m) 
and \>p_) — —e~ lip / 2 \fh') . On the contrary, for W < 
and \W\ > A one has \tp + ) = e~ iv/2 \rh') and |^_) ~ 
e iip/2 |^ Exactly at the resonance, W — 0, the eigen- 
states |^±) are superpositions of |m) and |m') with equal 
weights, and signs depending on the phase <p. 

Eq. (|45|l is sufficient to calculate the matrix element 



Vw 2 



(47) 



This result does not depend on the detailed form of Ha, 
provided that the condition of strong uniaxial anisotropy 



is fulfulled, and it is universal in this sense. It is large 
near the resonance, W 2 < A 2 , but it becomes small far 
from the resonance. This the only matrix element that 
plays a role in the relaxation for H = 0, as will be shown 
in Sec. IIV CI Also for H ^ in the case where the 
tunnel splitting A is due to the transverse anisotropy, 
this matrix element is dominant. In another particular 



case when 



H A ,S Z 



■ and thus the only source of 

the splitting A is the transverse field one has to take into 
account matrix elements \S± \ "0+) • This will be done 
in the next subsection. 



B. Beyond the two-state model: IM>_ \S±\ip + ) 
Consider the spin Hamiltonian Hs of Eq. I|15l) , where 



Ha satisfies 



Ha, S 2 



= and is strong in comparison to 



Hz- Due to this rotational invariance the orientation of 
the transverse field is unimportant, so we direct it along 
the x axis for simplicity: 



Hz — —g^B-HzSz ~ g^B-^xS x 



(48) 



and asssume H x > 0. For strong uniaxial anisotropy Ha 
one can obtain the splitting of resonance levels |m) and 
\m') {ml > to) perturbatively in H x : 



A 

where 



V, 



1 



m,m+l 



E m +i — E n 



-V, 



m+l,m+2 • 



Vrr, 



(49) 



V m , m +i = --^9Hb h x (m \S-\ m+ 1) . (50) 

In the particular case Ha = —DS 2 the calculation in Eq. 
(PI) yields^ 



A 



2D 



[(to' — to — 1)!] 



i(S + m')\(S-my. (g\i B H x 



(S-m')\(S + m)\ V 2D 



(51) 



One can see from this calculation that in our case in Eq. 
(Pl> is 



«p _ 



(52) 



for H x > 0. 

To compute matrix elements of S± between the tunnel- 
split states \ip±) one has to go beyond the two-state 
model of the preceding section. Taking into account that 
the states \rh) and \rh') are not pure |m) and |to') states, 
see Eq. (|41ll . one obtains small values for (rh\S-\m') 
and (rh' \ S+ \ rh) that are, however, essential in the spin- 
phonon relaxation in the case when A is caused solely by 
the transverse field. These small terms can be calculated 
perturbatively by building chains of elementary matrix 
elements that join the pure states, \m) and \m'). These 
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chains contain terms of H$ that do not commute with 
S z and the corresponding energy denominators, similar 
to the perturbative chain of Eq. (|49|l . The difference is 
that the elementary matrix element with the "external" 
operator S— does not contain — (l/2)gfi B H x , unlike all 
other elementary matrix elements, and the external oper- 
ator S- can be inserted into the chain atm'-m different 
positions. For so defined indirect (real!) matrix elements 
one obtains 



{m\S-\m') = (m'\S+\m) = (-1) 



m'-m+i ( m> - m) A 



9Vb h x 



(53) 

Note that this result does not use the specific form of 
Ha and it is thus universal, similarly to Eq. 1)47(1 . For 
S = 1/2 one has A = gn B H x and (S\S-\S) = 1, 
which is a correct result. Calculating the matrix elements 
between the tunnel-split states \ip±) given by Eq. 145(1 
with m => rh, m' => m', and ip = one obtains the 
unchanged result for (il>_ |5 , z | - + ) that is given by Eq. 
((47(1 . For transverse operators one obtains 



<VUS±IV>+) 



(m! — to) A 



±1 



W 



VW 2 + A 2 



Below we will need 



(rJ;_\S x \iP + ) = 



(m! — to) A 



W 



and 



(ip- \Sy\4>+) = ~i 



2gi-i B H x v'w 2 + A 2 
(to' — to) A 



(54) 
(55) 

(56) 



2gfi B H x 

Note that these results are only valid if A is due to the 



transverse field. When 



Ha, Sz 



^ the dominating 



source of A is Ha- In this case the matrix elements 
(ip_ \S±\tp + ^ are much smaller than (ip_ \S z \ip+) and 
they can be safely neglected. 



this property, together with the antiunitary character of 
K, leads to the relation 



(57) 



for the spin matrix element in the case of integer S. One 
can see that spin-phonon transitions are not ruled out 
by the time-reversal symmetry if the matrix element is 
imaginary. Complex conjugation in Eq. 1(57(1 makes the 
situation different from the cases of spatial symmetries. 

A specific example showing the absence of a time- 
reversal selection rule is the crystal field Hamiltonian 
Ha from Eq. J2J). The contribution to H s -ph from the 
phonons rotating the lattice around z-axis is imaginary 
and is given by Eq. ((B5(l . On the other hand, eigenfunc- 
tions \ip±) of Eq. 1(45(1 are real for E > [see comment 



below Eq. |B4|)]. Thus (ip_ 



and Eq. ((571) is satisfied by (-0 



P h 



-ph 



is imaginary 



In the lattice frame, as well as in the laboratory frame, 
if our new method is used, one has to calculate matrix 
elements of the spin operator S [see Eqs. (127}) and J3§J] 
between the states \ip + ) and \ip_). I n contrast to the 
spin part of iJ s _ p h, the operator S breaks time- reversal 
symmetry, see Eq. I(C3(I . Thus one obtains the relation 



(</UShA + }HV>_|ShA + >* (58) 

instead of Eq. ((57(1 . For the biaxial model with E > 
the states and \ip__) are real, and one obtains the 

selection rule (ip_\S y \ip + ) — because S y = (i/2)(S- — 
S+) is imaginary. On the other hand, S x and S z are real, 
and time-reversal symmetry does not lead to selection 
rules for them. However one obtains (ij}_\S x \xf! + ) = 0, 
too, that can be shown with the help of S x — (1/2) (S- + 
S + ) and Eq. J25J|. Thus the matrix element (ip_\S x \i/} + ) 
given by Eq. I(47l) , corresponding to phonons rotating the 
lattice about the z-axis, is the only matrix element that 
is responsible for the relaxation between the tunnel-split 
states of the spin Hamiltonian for H = 0. 



C. Role of time-reversal symmetry 

Let us now discuss the role of time reversal symmetry 
for spin-phonon transitions. For the total Hamiltonian H 
to be invariant under time reversal, the external field H 
must be zero. Thus we will consider only this case here. 
For H = 0, tunneling can only arise from the transverse 
anisotropy which lifts the degeneracy of the eigenstates 
| to) and | — to) of the longitudinal part of the crystal 
field Hamiltonian Ha- According to the Kramers' theo- 
rem, this degeneragy is lifted only for integer spins S. Let 
\ip±) denote the corresponding tunnel-split eigenstates of 
Ha- It will be shown in Appendix IClthat \ip±) are eigen- 
states of the time reversal operator K with eigenvalues 
±1, i.e. \ip + ) and have opposite parity with respect 
to time reversal. The spin-lattice Hamiltonian iJ s _ p h in 
the laboratory frame that is given by Eq. ((12(1 is invariant 
under time reversal. It will be proven in Appendix IClthat 



V. SPIN-PHONON RELAXATION FOR 
TUNNEL-SPLIT STATES 

The rates of spin-phonon transitions can be calculated 
with the help of the Fermi golden rule. One should distin- 
guish between two situations: when the tunnel splitting 
A is caused by the terms in H a that do not commute 
with S z , such as the transverse anisotropy in Eq. J2J, 
and when A is caused by the transverse field in Hz- The 
physical difference between these two cases is that Ha is 
defined in the lattice frame and thus it is rotated by the 
transverse phonons, whereas Hz is defined in the labora- 
tory frame and it is not rotated. We will see that in the 
first case, in the absence of the field, relaxation is due 
to the phonons rotating the lattice around the z axis. 
To the contrary, in the second case (Ha commutes with 
S z ) these phonons are decoupled from the spin and they 
produce no effect. 
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A. Tunneling induced by the anisotropy 

In this case, as was shown at the end of Sec. lIVCl the 
transition matrix element in Eq. I|36[l is due to the oper- 
ator S z only for H = 0. For nonzero fields, S z provides 
the dominant contribution to the matrix element because 
\S Z \ 4>+) ~ 1 in the vicinity of the resonance [see Eq. 
(|47fl ] , whereas other matrix elements can be shown to be 
small. Thus one can write 



s— ph 



(59) 



With the help of Eq. JTTJ and Eq. J2gJ, Eq. (JSDJl can be 

rewritten as 



\H R . 



ph 



E 

kA 



-kA 



«k 



(60) 



where we used the designations of Eq. (|25[1 and 



14a 



V8Ml^kA 
«3Mb [H x [k x e kA ]]J 



{fuv [k x e kA ]. 



(61) 



The decay rate W! |_ of the upper spin state into the 

lower state, accompanied by the emission of a phonon, 
and the rate W-\ of the inverse process are given by 



where n u 



W+- 



W 



+ 1 



(62) 



1) with = l/{k B T) is the 



phonon occupation number at equilibrium and 
Wo = -^^|^kA| 2 27n5(^ kA -cjo) 



(63) 



kA 



The master equation for the populations of the spin states 
n+ and n_ satisfying n+ + n_ = 1 and 



n + = -W- + n + + W + -n- 
= -Tn+ + W+- 
defines the relaxation rate 



r = w+- + W- 



W (2n u 



1) 



(64) 



(65) 



One can see that the two terms Eq. ()61|l do not inter- 
fere and that only transverse phonons, A = t, are active 
in the relaxation process. In the first term of Eq. I|61|) 
one can use [k x e ki ] = ±&e k t', where t and t' denote 
different transverse phonons. Summation over polariza- 
tion vectors of transverse phonons and averaging over the 
directions of k can be performed with the help of the for- 
mulas given in the Appendix iDl Replacing TV -1 ^ k . . . 
by vq J d 3 fc/(27r) 3 . . . (vq being the unit cell volume) and 
using Eq. (|47(1 one obtains 



W 



[(to'-to)/2] 2 A 2 ^ (^ 2 + 7 2 g 2 ) 

12-Kh Mv\ lu 3 d 
[(to' - to)/2] 2 A 2 uj (luI+^HI) 



12-kK 



(66) 



1 /3 

where p is the mass density, u>d = v t/v is the Debye 
frequency for the tranverse phonons, ujq is given by Eq. 
PH>. and H\ = H 2 + H 2 . In Eq. ® the term w 2 , in 
the brackets is due to the transverse phonons that rotate 
the lattice around the z axis, whereas the transverse-field 
term is due to the phonons that rotate the lattice around 
x and y axes. 

The beauty of Eq. I|66|) is that it gives a universal ex- 
pression for the transition rate, which does not depend 
on the exact form of the crystal-field Hamiltonian, pro- 
vided that the uniaxial anisotropy dominates, and is ex- 
pressed entirely in terms of independently measurable pa- 
rameters. Equations (|39fl . I|44|) . and l|66|l show that for 
(to' — to) 2 ^> 1 the contribution of the longitudinal bias 
field to the relaxation is much stronger than the contri- 
bution of the transverse field. 



B. Tunneling induced by the transverse field 

The spin matrix elements for this case have been com- 
puted in Sec. liVBl The calculation of ft |_ is more com- 
plicated than in the previous subsection because one has 
to take into account the transverse components in the 
field term in Eq. , 

0-+ = {iuoe x + jH z e y ) 54>_ + x 

+ (iLu e y - jH z e x + -fH x e z ) S(j)_ + y 

+ (iuj e z - yH x e y ) 84>_ + z . (67) 

The transition matrix element, Eqs. Ij27(l or Ij36(l . can be 

presented in the form 



P h 



K 



x"Y — u x 



-K 



(68) 



Using Eqs. I|47(l and (|56|l one obtains 

K z = -i (ip_ \S z \ip + ) tihj + (ip_ \S y \ ip + )gfi B H x = 0. 

(69) 

This result was expected because phonons rotating the 
lattice around the z axis cannot cause any effect for 

0. To the contrary, K x and K y are nonzero: 



Ha, S z 
K x = 



-i (ip_ \S x \ip + ) huj Q - gp, B H z (ip_ \S y \ip + ) 
i (to' — to) A 



2gp B H x 



(W-gp B H z ) 



(70) 



and 



K y = -i(ip_\Sy\ip + )tuJo + gp, B H z (ip_\S x \tp + ) 
-gp B H x (ip_ \S z \ip + ) 
, (to' — to) A 



-gp B H x S- 



-A 



| W{W~gp B H z ) 

{g^ B H x f 



where A is given by 

A = 1 



(gp B H x ) 



(71) 



(72) 
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and W is related to H z by Eq. 

Quantization of the phonon field, Eq. i|ll|) . yields Eq. 
(EOJ, where 



Vi 



kA 



ji^ [k x e kx ] x + K y [k x e kA ] y | 



(73) 

The relaxation rates between the states m and to' are 
given by Eq. H62|) . where 



Wo 



1 iKr 



127l7l 



Mv 2 



(74) 



and X x and K y are given by Eqs. I|7U|I and (|7f |l . respec- 
tively. The general result is cumbersome but it simplifies 
for the ground-state resonance, m — —S and to' = S: 



W 



S 2 A 2 u (~fH ± y 
12nhMvf uj 3 d 



(75) 



where we have replaced H x =>■ H± , the field perpendicu- 
lar to the anisotropy axes, and 



Q 



l - 



/ ^0 



1 / W 



2S \gn B Hj 



1 - 



(- 



W 



2S J \g^ B H A 



(76) 



In the case of S — 1/2 one has A = g/.i B H±, so that 
Q and Wo are zero. This is expected as the crystal-field 
anisotropy Ha disappears for 5 = 1/2 and phonons do 
not couple to the spin. In general, for the ground-state 
splitting one has A oc H^ s and A <C g/J, B H± for any 
S > 1/2. 



VI. DISCUSSION 

We have studied a universal mechanism of the relax- 
ation of a spin state in a solid in the presence of the 
magnetic field. It corresponds to the generation of the 
elastic twist mandated by conservation of energy and to- 
tal angular momentum. This mode of relaxation must 
be the dominant one for a rigid spin cluster embedded 
in an elastic medium. Simple universal expression, Eqs. 
(|27|l and (|36|l with Eq. (|28|l , for the transition matrix ele- 
ment has been derived for an arbitrary spin Hamiltonian 
and an arbitrarily directed magnetic field. The method 
of computing the transition matrix elements, presented 
in this paper, consists of two steps. In the first (needed 
by any theory) step one must obtain the eigenstates of 
the spin Hamiltonian, the transition between which is go- 
ing to be studied. Once these eigenstates are known, the 
computation of the spin-phonon transition rate by our 
method reduces to the calculation of the matrix elements 
of S between the eigenstates of interest. We have studied 
in detail the transitions between tunnel-split spin states, 
which arc difficult to compute by conventional methods. 
The role of the time-reversal symmetry has been ana- 
lyzed. As demonstrated in Sec. II V CI not all transitions 



between spin states of different parity with respect to 
time reversal are ruled out by time-reversal symmetry. 
This is at variance with the case of spatial symmetries. 

Universal formulas have been obtained for two partic- 
ular cases of the strong uniaxial anisotropy. The first 
case is that of the tunnel splitting created by the crys- 
tal field, Eq. I|66|) . The second case if that of the tunnel 
splitting created by the transverse magnetic field in the 
presence of an uniaxial magnetic anisotropy, Eq. i|74[l . 
The first case corresponds to, e.g., the intensively studied 
S = 10 Fes molecular cluster, while the second case cor- 
responds to, e.g., the recently studied S — 4 N14 molec- 
ular cluster^ For both cases, we have obtained relax- 
ation rates, and their field dependence, in terms of inde- 
pendently measurable constants and with no adjustable 
parameters. Among other applications, these rates are 
responsible for the linewidths measured in electron spin 
resonance experiments. Our results opens the way for 
an accurate comparison between theory and experiment 
on the field dependence of the relaxation rate. Among 
many possible applications our theory also has an im- 
portant consequence for the industry of spin-based solid- 
state qubits. It provides parameter-free, mandated by 
symmetry, lower bound on the decoherence of any such 
qubit. 

While multiphonon processes can be studied by the 
same method, this paper is limited to the rates of direct 
one-phonon processes that dominate spin-lattice transi- 
tions at low temperature. Two-phonon Raman processes 
that contribute at higher temperatures will be studied 
elsewhere. 
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APPENDIX A: SPIN-PHONON RELAXATION 
FOR ADJACENT LEVELS 



Here we shall study the spin Hamiltonian 
H S = H A + H Z = -DS 2 - g[i B B z S z . 



(Al) 



The exact energy levels of this Hamiltonian are given 
by Eq. (|3"?|) with Em = —Dm 2 . We consider the spin- 
phonon relaxation between the adjacent levels of H$, m 
and to' — m + 1. The energy difference between these 
levels is 



W = hu = E m - E m+1 
= D(2m+ 1) +gfj, B H z 



(A2) 



(we assume W > 0). We first study the (to, to + 1) tran- 
sitions by a conventional method that employs tilting of 
the anisotropy axis by transverse phonons, see, e.g., Refs. 
IT6II24I In this method one writes Ha in the form 

H A (n) = -D(n-S) 2 , (A3) 
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where n is the direction of the anisotropy axis, 

n = e z + Sn, Sn = S<j> x n. (A4) 

Expanding TJ^n) up to the linear terms in S<f) one ob- 
tains -ff^n) = —DS 2 Z + -ff s -ph with 

4-ph = -D [(Sn ■ S) (e,-S) + (e,-S) (Sn ■ S)] = 
-D [S X S Z + S Z S X ] 5^ y + D [S y S z + S z S y ] Sfa. (A5) 

The same result follows from the calculation of the com- 
mutator in Eq. i|12fl . We shall calculate the transition 



matrix element between the states yf! 
and |^_) = \m + 1) ® |n k + !)• One obtains 



+/ 



\m) |nk) 



s — ph 



*+) = — l m ,m+i (2m + 1) (n k +l|^_|n k ), 

(A6) 



where l m , m +i = y/S(S + 1) — to(to + 1) and 

<S0_ = - iS(j) y . (A7) 

Note that the contribution to Eq. I|A6(1 of the transverse 
phonons of definite chirality, Eq. I|A7(1 . having l z = 1 
projection of the angular momentum onto the z-axis, is 
in accordance with the conservation of the total angular 
momentum for the (to + 1, m)-transition. 

Quantization of lattice rotations with the help of Eq. 
ifTTJl yields Eq. (JUJ with [cf. Eq. 



kA 



m,m+l 



(2m + 1) 



■ikr 



k x e kA ] ■ (e x - ie y ) 
y/8MfkJ U x 



The relaxation rates between the states to and to + 1 are 
given by Eqs. ^ and ||BSJ with V4 A defined by Eq. (jUJ. 
With the help of the formulas listed in the Appendix iDl 
one obtains 



Wo = 



(2TO + l) 2 lm,m+l D 2 ^m,m+l 



2AttH 



Mv 2 



(A9) 



The limit of D — * corresponds to a free spin. In that 
limit, although the states \m) and |to + 1) are still sep- 
arated in energy due to the magnetic field, the rotation 
of the lattice cannot cause any relaxation and Eq. (|A9(I 
produces zero result. 

It is instructive to see how the above result can be 
obtained by our new method using Eqs. (|36l) and 

with f2 |_ = (n k + l|n|rik). As the transition \m) — > 

|to + 1) is due to the spin operator S+, one obtains 



s— ph 



4-, 



-lm,m+i(nk + l|0_|n k ), (A10) 



APPENDIX B: ROLE OF ROTATIONS AROUND 
THE z AXIS 

In Sec. IV Al we have calculated the spin-phonon relax- 
ation rate between the tunnel-split resonance states of 
the spin in the case when the tunnel splitting A is due 
to the terms in the crystal-field Hamiltonian Ha that do 
not commute with S z . The result, Eq. (|66|) . is universal 
and it does not depend on the detailed form of Ha- We 
have seen that the relaxation is generated by the phonons 
rotating the lattice around the z axis. The physics of this 
process is transparent: in the rotating frame the spin feels 
the magnetic field along the Z-axis. This field couples to 
S z , which produces the relaxation. In this Appendix we 
will obtain this result for a biaxial spin model by the 
conventional method used in the previous Appendix. We 
shall demonstrate that in the conventional method the 
relaxation arises from the rotation of the hard axis by 
the phonons. As we shall see, the computation of this 
effect, even for the simplest spin Hamiltonian, is signif- 
icantly more cumbersome that the computation of the 
transition rate by our new method. 

Consider the spin Hamiltonian H$ of Eq. I|15l) with Hz 
given by Eq. J^J) and 

H A = -DS 2 Z + E (S 2 X - Si) = ~DS 2 + \e (S\ + S 2 _) . 

(Bl) 

We assume E > for simplicity. If E is negative, one 
can change its sign by rotating the coordinate system 
by 7r/2 around the z axis. For E <C D the energy lev- 
els of this system can be approximately characterized by 
the quantum number to. Any two levels to' > to can be 
brought into resonance by the longitudinal component 
of the magtetic field H z , see Sec. IIVI The tunnel split- 
ting A can be calculated with the help of the high-order 
perturbation theory: 



-V, 



m-f2,m+4 ■ 



V, 



m' — l.ra' 



A -v 1 

„ — Vm,m+2 „ j-, 

where 



V m . m+ 2 — -E (to \S 2 _ J to + 2) — -El m ^ m+ il m+ i 



(B2) 



m+2- 

(B3) 



Calculation in Eq. I|B2|I yield; 

2D 



[(m , -m-2)\\y 



'(S + to')! (S -to)! 



(S 



(s- 



2D J 



(m — m.) /2 



.(B4) 



where 57 _ 



iCl y . For H = H z e z this gives 



pii 



2^m : m+l (fi^o ~ 9^B H z) («k + 1 1<5</>_ |n k ) . (All) 



Here the magnetic field cancels according to Eq. 
and one obtains Eq. (|A6p . 



One can see from this calculation that in our case in Eq. 



(Pl> is 



1. Let us now consider the phonons that 



rotate the lattice around the z axis. It is easy to obtain 
that the result of such rotation is 



-ph 



-iE(S\ - S 2 _) 54> z = VS<t> z 



(B5) 



Now, similarly to the calculation in Sec. IIVBI one can 
compute the matrix element of _ff s _ p h between slightly 
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delocalized states \fh) and \m') [see Eq. (I41H ] by building 
a minimal perturbative chain between the pure states |m) 



and 



and identifying the result with A given by Eq. 



(|B2ll . As Hs-ph can be inserted at (to/ 
in the chain, one obtains 



to) /2 positions 



(in \S 
and thus 



2 'to' 



) = (m'\S\\rh) 



A 



E 



(B6) 



(in 
(in 1 



m — m . 
to' — m 



-A. 



(B7) 



For the spin eigenfunctions given by Eq. I|45[) with e %v 
lone obtains 



v 



The spin-phonon matrix element is then given by 



- P h 



(B8) 



(B9) 



The latter coincides with the result obtained by our new 
method, Eq. (SgJ with \S Z \ given by Eq. I|47|) 

and Ml |_ jZ = iMjQ(f> , z . The advantage of the new 

method is apparent as it only involves the trivial compu- 
tation of the matrix element of S z . 

To gain a deeper insight into the difference between the 
two methods, one can start with the general form of the 
eigenstates \^±) given by Eq. (|4*U|l and compare general 
expressions for the spin-phonon matrix element due to 
the rotation around the z axis, used by the two methods. 
The traditional method uses the worked out form of the 
commutator in Eq. I112H while the new method uses Eq. 
( 1351) . The equivalence of the two results relies on the 
identity 



m+2 



A 
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s 



C*mC. 



(BIO) 

where A = E + — EL and the matrix element V m ^ m +2 is 
due to the terms in Ha that do not commute with S z . 
(We consider for illustration the model with V m . m +i = 
and zero bias, H z = 0.) The left-hand side of Eq. 
(|B10p corresponds to the worked out commutator and 
it vanishes if one approximates \ i f?±) by the two-state 
model of Sec. IIV Al i.e., if one neglects all C± >m i' with 
m" =/= m,m'. Leaving C± m with all m leads to the 
correct result of order A in the left-hand side of Eq. 
(|B10p . but only after one accurately accounts for the 
cancellation of many small terms containing powers 
of E/D <C 1. On the contrary, in the right-hand side 
of Eq. (|B10|) the small quantity A has been already 
factored out, and it is sufficient to make the two-state 
approximation in the sum over to. 



APPENDIX C: ROLE OF THE TIME-REVERSAL 
SYMMETRY 

In quantum mechanics time reversal is represented by 
the operator K that satisfies K~ x = K\ and, in addi- 
tion, performs complex congugation^i The latter comes 
from the requirement that the lhs, ihdt \ip{t)} of the 
time-dependent Schrodinger equation is covariant under 
t — > — t (which changes i — > —i), only, if one takes its 
complex conjugate. The latter makes K antilinear: 



K [ci | Vi) + c 2 | v 2 >] = c\K IVx) + c* 2 K |v 2 



(CI) 



Because of the antilinearity, one has to specify whether 
K acts to the right (by default) or to the left: 



K 



1> 



(jr|p)) = ((^|£)|v)*- (C2) 



Operator K is called antiunitary because of its antilin- 
earity and its property K^ 1 = . 

Transformation of the spin operator S under time re- 
versal follows from the requirement that S should behave 
as the orbital momentum L and it thus should change its 
sign under time reversal: 



s' = ksk f = ksk- 1 



(C3) 



An arbitrary spin state \ip) can be represented as a su- 
perposition of the eigenstates \m) of S z . Choosing the 
phase factor of \m) to be one, i.e., taking |m) as real 
makes the matrix elements (m\S z \ m!) and (to \ S± \ to') 
real. In that case the explicit form of the time-reversal 
operator is K — QU = UQ, where U = exp (— inS y ) is 
a unitary operator and Q makes complex conjugation. 
Indeed, U changes the signs of S x and S z , whereas Q 
changes the sign of S y , so that Eq. (|C3I) is fulfilled. 

In order to determine the action of K on the tunnel- 
split spin eigenstates \ip±) of the crystal- field Hamilto- 

nian Ha we have at first to determine the action of K on 
\m). Application of KS Z = -S Z K and KS+ = -S-K 
[that follow from Eqs. I|C3|) and IClfl ] onto |m) yields, up 
to an irrelevant global factor, 



K\m) 



-1) 



S-mi 



(C4) 



As a consequence one obtains K 2 \m) = (— 1) 2S |to) , i.e., 
K 2 = 1 for integer spins S. 

In zero field H$ reduces to Ha that must be time- 



K,H A 



0. 



reversible, KH A K 1 = Ha, i.e., it satisfies 

For integer S the tunnel-split eigenstates \il>±) of Ha are 
nondegenerate, thus they must be eigenstates of K: 



K \ip ± ) = e± \tp±) 



(C5) 



with eigenvalues e±. Since K 2 = 1 for integer spin S", it 
must be e± € { — 1,1}. Let us prove now that £ + and 
e_ have different signs. To this end, we decompose Ha 
into the longitudinal and transverse parts, H? ng and 
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^trans 



long 



and thus 



The former satisfies S z , H 
it has | ± m) as at least twofold degenerate eigenstates. 
The transverse part H A rans has matrix elements between 
different |m) and thus it removes the degeneracy between 
| ± m) for integer S. Let us introduce 



H A (\) 



= H l ° ns 



^j^- trans 



(C6) 



with A real. In the limit A — > one can find V±) analyt- 
ically. Using Eqs. (|C4|) and (|C1|I . as well as the relation 
(— l) 2m = 1 for integer S, one can check that 



|V±(0)) = -p [e* x |m) ± e~ lx | - m)] 
V2 



satisfy Eq. (|C5|) with 

e±(0)=±(-l) s - m . 



(C7) 



(C8) 



Since //a (A) is continuous in A, the eigenstates |^ ± (A)) 
and, in turn e±(A), must be continuous, too. This con- 
tinuity and the discreteness of e± implies that e± are 
independent of A. Thus one obtains £+(A) = — £_(A) for 
all A, including A = 1. Thus \^p±) = |^±(1)) inherits its 
parity from the unperturbed eigenstates. Consequently, 
we get an interesting result that the tunnel split eigen- 
states of an arbitrary crystal field Hamiltonian must have 
opposite parity with respect to time reversal symmetry, 



(C9) 



The spin-phonon Hamiltonian i? s _ p h is invariant un- 
der time reversal, i.e., H!,_ ph = KH s - p hK^ = H s - P h- 
One can obtain a symmetry relation for the spin matrix 
element if one makes time reversal of both ff s _ p h and the 
spin states. Inserting K*K = K~ l K = 1 into the matrix 
element one proceeds as follows: 



1p_ i? s _ph 1p+ 



4>- 



s— ph 



K 



e+e- (ip 



H, 



s— ph 



H. 



s— ph 



(CIO) 



where we have used Eqs. l|H2jl . (tC5l . and JHH). This 
proves Eq. I|57l) . Spin operator changes sign under time 
reversal, Eq. (|C3|) . thus a similar procedure leads to Eq. 

APPENDIX D: USEFUL RELATIONS 

To calculate the spin-phonon relaxation rates, the fol- 
lowing useful relations can be used. 

Summation over the two transverse polarizations: 



[k x e ktl ] = ±ke 



kt 2 



(Dl) 



and 



E(k ■ a) (k ■ b) ,„ , 
(e kt ■ a) (e kt ■ b) = (a - b) - * ^ '-. (D2) 



i=ii,t 2 



Averaging over the directions of the vector k: 



((k-a) (k.b)) = T (a-b) 



(D3) 
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